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Conventional electronics are based invariably on the intrinsic degrees of freedom of an electron, 
namely, its charge and spin. The exploration of novel electronic degrees of freedom has impor- 
tant implications in both basic quantum physics and advanced information technology. Valley as 
a new electronic degree of freedom has received considerable attention in recent years. In this 
paper, we develop the theory of spin and valley physics of an antiferromagnetic honeycomb lat- 
tice. We show that by coupling the valley degree of freedom to antiferromagnetic order, there is 
an emergent electronic degree of freedom characterized by the product of spin and valley indices, 
which leads to spin-valley dependent optical selection rule and Berry curvature- induced topologi- 
cal quantum transport. These properties will enable optical polarization in the spin- valley space, 
and electrical detection/manipulation through the induced spin, valley and charge fluxes. The 
domain walls of an antiferromagnetic honeycomb lattice harbors valley-protected edge states that 
support spin-dependent transport. Finally, we employ first principles calculations to show that the 
proposed optoelectronic properties can be realized in antiferromagnetic manganese chalcogenophos- 
phates (MnPXs, X = S, Se) in monolayer form. 



The exploration of novel electronic degrees of freedom 
[T]-[4] has been a fairly important topic recently, as the 
future of information technology is likely to be levered 
upon it. The intrinsic degrees of freedom of an elec- 
tron, namely, its charge and spin, has been the basis for 
the society-transforming information technologies, i.e., 
electronics and spintronics. Additional electronic degree 
of freedom, if present, will offer immense potential for 
information encoding and manipulation at the micro- 
scopic level. The notion of valleytronics [1 on honey- 
comb lattices has received considerable attention in re- 
cent years. [SHTO] When the centrosymmetry of the hon- 
eycomb lattice is broken, such as in gapped graphene, 
there arise an inequivalent, degenerate pair of valleys in 
the momentum-space electronic structure. The valley 
excitations, protected by the suppression of intervalley 
scattering, have contrasting optical and transport prop- 
erties ensured by quantal helicity. O [6] In the most recent 
experimental progress in monolayer group VI transition 
metal dichalcogenides, the identity of valleys manifests as 
valley- selective circular dichroism, leading to substantial 
valley polarization with circularly polarized light, offer- 
ing a potential arena to the eventual realization of val- 
leytronics. [ilElIIo] 

Central to these endeavors are two tasks, namely, de- 
veloping theoretical paradigms and subsequent materi- 
als discovery, the latter of which makes possible exper- 
imental measurements that put the theory to test 
[3[9l[T0]. Evidently, it is important to broaden the choice 
of materials, beyond transition metal dichalgogenides, [7^- 
[TO] with which novel degrees of freedom (beyond charge 
and spin) of Bloch electrons can be accessed. In previous 
theoretical and experimental developments, attention has 
been paid to the absence of inversion center in the lattice 
space group. Here we show that the pseudospin symme- 
try present in the initially symmetric honeycomb lattice 



hints at a nontrivial transformation, which leads to an 
emergent degree of freedom characterized by the prod- 
uct of spin and valley indices. We further propose that 
such a Fermionic system may be observed on a bipartite 
honeycomb lattice that assumes the Neel antiferromag- 
netism (afm), which can possess chiral electronic exci- 
tation concomitant to the spin density wave, as well as 
spin-dependent transport properties on the edge states. 
In the last part of the paper, we employ first principles 
calculations to show that the proposed physics can be 
realized in antiferromagnetic manganese chalcogenophos- 
phates (MnPXa, X = S, Se) in monolayer form. We also 
analyze the general consequence of symmetry in the afm 
Fermionic honeycomb lattice. A potential transition to 
the topologically non-trivial quantum spin Hall state is 
briefly discussed. 

Valley as an electronic degree of freedom has been 
suggested in a spinless Fermion model on a honeycomb 
lattice with a symmetry-breaking perturbation to the 
sublattices.|5l E] A honeycomb lattice can be defined as 
a three-connected two-dimensional net, with the connec- 
tion vectors for bond length h pointing toward three near- 
est neighbors, di 2 = (±a/3x + y)6/2, and da = —by. 
The corresponding lattice constant a = VSb. The two 
sublattices of a honeycomb structure is conventionally 
denoted A and B, respectively, which correspond to a bi- 
nary degree of freedom called isospin. The low-energy 
quasiparticle states at K± = ±47r/3ax are assigned a 
valley index r = ±1, indicating the valley contrasting 
physics. It is tied to physical measurables, such as the or- 
bital magnetic moment, A4{'K±) = r/ig, where /ig is the 
effective Bohr magneton. [5 Concomitant to the orbital 
magnetic moment is non- vanishing Berry curvatures, pJJ 
^{K±) = hM{K±)/ee{K±), where e and e{K±) are 
the electronic charge and the band energy at K± , 
respectively. ^ 
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When the spin degree of freedom {s = ±1/2) is taken 
into account, spontaneous symmetry breaking becomes 
viable in the dichromatic Shubnikov group. Imposing 
dichromatic coloring on the spin and lattice, we propose 
a spin-full Hamiltonian for the low-energy quasiparticles 
near K± of a honeycomb lattice, 

H^^^^ = VYSoir^a^Px + TQCFyPy) + ms^rocr^, (1) 

where p is the momentum and is the massless Fermi 
velocity. Here, r^^ (J a and (o^ = ^, ^, ^, 0) are the Pauli 
matrices for the valley, isospin and spin degrees of free- 
dom, respectively. The mass term, m, corresponding to 
a symmetry-breaking perturbation, admits a band gap 
A = 2m for both spins. With the spin-dependent Hamil- 
tonian, M(K±) = 25r/ig. Therefore, the product of spin 
and valley indices identifies a new degree of freedom of 
electrons, which we shall call the coupled spin- valley, s-r, 
degree of freedom. The s • r index bears the virtue of a 
good quantum number, for so long as inter- valley scat- 
tering is suppressed it mimics the two spin states of an 
electron. 




FIG. 1: The antiferromagnetic honeycomb lattice, (a) Neel 
antiferromagnetism of a honeycomb lattice. Up and down ar- 
rows indicate a staggered spin-density wave, (b) A schematic 
of the spin-dependent lattice potential corresponding to the 
afm order, (c) Low-energy quasiparticle bands of a NN hop- 
ping afm Hamiltonian, Solid lines assume zero spin- valley cou- 
pling, while dashed lines take into account the spin- valley cou- 
pling. Inset: the Brillouin zone and high symmetry points. 

As the spinless Fermion model describes a non- 
magnetic system, [5 the dichromatic Hamiltonian of 
equation ([T]) represents an antiferromagnetic system. In- 
spection of Hamiltonian ([T]) reveals that the s-r degree of 
freedom arises upon engaging the electrons with the spin- 
dependent on-site Hubbard potential, Un^^^^n^- ^\ which 
characterizes the Coulomb interaction of opposite spins 



on site j = A, B. [12] Within the Hartree-Fock approxi- 
mation, the on-site Coulomb energy of two spins can be 
decoupled as m^^^n^-*^ +m^~^^n^- ^\ upon an immaterial 
energy shift. We define a spin-dependent mass term, 

m^'^ = 2sm, (2) 

where 

^=j\{n^r^)-{4^)l (3) 

characteristic of an afm spin-density wave (Fig. la). It 
is now apparent that with the afm order, each spin sees 
different potentials on the two sublattices, as represented 
by m*^*^ (Fig. lb) and described by equation ([T]). 

To pinpoint the essential physics of spin and val- 
ley on an afm honeycomb lattice, we analyze a tight- 
binding Hamiltonian invoking the spin-dependent mass 
term above and a nearest-neighbor (NN) hopping, t, 

= t'l^sodx + t'^sody + mSzCTz, (4) 

which in the neighborhood of K± can be linearized to 
equation ([T]). The NN hopping is given as 

3 

= ^ exp(-zk • d^) = tk + i^k, 

where are the vectors pointing toward the three NN's. 
The band structure is shown in Fig. Ic. The bands are 
spin degenerate at each k-point. Band gaps A = 2m 
are indeed opened at the valleys. An ad hoc spin- valley 
coupling may be introduced with the parameter J, as 
H^^^ = dszTzCFz- The spin- valley coupling preserves the 
spin degeneracy, but leads to a renormalization of the 
valley gaps. The band gaps become enlarged at one val- 
ley and reduced at the other, that is, A = 2{m — T5). The 
fact that spins remain degenerate in this system has to do 
with the invariance of the Hamiltonian under simultane- 
ous time reversal (T) and spatial inversion {V)^ although 
neither T nor V alone commutes with the Hamiltonian. 
We shall revisit this symmetry, O = VT^ later in the 
paper. 

A few interesting experiments become immediately 
compelling, to probe and manipulate the s • r degree 
of freedom. When the spin-valley coupling is absent or 
weak, the two valleys can be considered degenerate (Fig. 
2a). The system will have s • r-selective circular dichro- 
ism (CD), similar to gapped graphene[6i and monolayer 
M0S2 [7WT0 . However, because the two spins adopt op- 
posite masses, the optical selection rule at becomes 
7^(5, r) = 25r, where r] is the degree of circular polariza- 
tion mil]. For example, when left-polarized light (cr+) 
illuminates the sample, spin up electrons are excited to 
the conduction band at K+, while spin-down electrons 
at K_ , as illustrated in Fig. 2a. When a right-polarized 
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FIG. 2: Spin and valley physics of an afm honeycomb lattice, 
(a) s • r-selective circular dichroism, in the absence of spin- 
orbit coupling, (b) When spin-valley coupling is present the 
valleys can be doped asymmetrically, (c) Electron spin (bot- 
tom of conduction bands) fluxes under the action of Berry 
curvature of the Bloch bands and in-plane electric field. The 
spin-up and spin-down currents show in blue and green, re- 
spectively. Solid and dashed lines stand for the currents from 
K+ and K_, respectively. The spin and valley, (s,t), indices 
are indicated in the parentheses. E is an applied in-plane 
electric field. 



photon is absorbed, it excites down-spins of and an 
up spins of K_. 

Berry curvature also induces topological quantum 
transport that allows for electric detection of 5 • r po- 
larization, via the anomalous velocity of Bloch electrons, 
E X l^k. [11] In the above dynamical polarization 
of charge carriers, the spin or spin holes under the ac- 
tion of Berry curvature, l](k, 5), of the Bloch bands will 
exhibit circular dichroic Hall effect (CDHE). Because of 
the relation with orbital magnetic moments, both Berry 
curvature and transversal conductivity {(Jxy) depend on 
s • r. Charge carriers with s-r = ±1/2 (excitable by left- 
and right-polarized lights, respectively) will have oppo- 
site transversal conductivity, as shown in Fig. 2c. There- 
fore, CDHE is a non-equilibrium charge Hall effect in the 
presence of a circularly polarized radiation field. 

If we dope the system with electrons or holes at equi- 
librium, then an applied in-plane electric field will drive 
a transversal "valley" current. Spin currents from all 
s-r = ±1/2 contribute to the transversal transport, 
resulting in a net accumulation of valley moments and 
orbital magnetization at the upper and lower edges of 
the sample (Fig. 2c) with zero transverse charge current. 
This is the valley Hall effect. In the case of strong spin- 
valley coupling, the gaps at K± become different, and 
thus have different levels of doping. Suppose we dope 
the system, say, at with electrons (Fig. 2b). Under 
an in-plane electric field, the carriers will produce a net 
transversal spin current without charge current, that is. 



the pair of currents with (5,r) = (±1/2,1) in Fig. 2c. 
This is the anomalous spin Hall effect (SHE). When we 
switch from n-doping to p-doping (Fig. 2b), the SHE 
will be characterized by a transversal spin current in the 
opposite direction, arising from spin holes. 




FIG. 3: Spin-polarized edge states, (a) The band structure 
in the presence of a zigzag magnetic domain on the afm hon- 
eycomb lattice, derived from NN hopping Hamiltonian. A 
model of a zigzag domain wall is used, which extends 50 unit 
cells on both sides. The bulk states are lumped into shaded 
blobs. The mass is set to m — t/2. The edge states of two 
spins are shown in blue and green, (b) The wavefunctions of 
two spins at the same Dirac points. The amplitudes are con- 
volved on 2-dimensional Gaussians centered on lattice sites 
for visualization (scale bar in arbitrary units). 

Also to be assessed is the topological domain walls 
that occur naturally with the afm order, of which the 
band structure is shown in Fig. 3a. There are two spin- 
polarized bands that arise at the domain wall and inter- 
sect to yield a pair of Dirac points, while the bulk states 
are always spin degenerate. The spin and momentum 
are clearly locked at individual Dirac point. Unlike the 
edge states in a topological insulator, these edges states 
cannot offer the spin-selective channels. However, if a 
spin is injected into the edge with a prescribed direction 
of momentum, it can migrate ballistically, enjoying the 
suppression of back scattering offered by valley protec- 
tion. It is also of interest to note that the wavefunctions 
(Fig. 3b) at the same Dirac points have opposite parity 
for opposite spins. The odd and even parity boundary 
state may be utilized as spin and/or valley filter, as linear 
defect in graphene. [13] 

Taking one step from the above theoretical exposition, 
here we suggest actual materials in which the proposed 
spin and valley physics can be observed. The selection 
of materials should meet a few criteria. The candidate 
materials must 1) have characteristic afm order on a hon- 
eycomb lattice; and 2) be a semiconductor, with direct 
band gaps at high-symmetry K+, K_; and 3) have in- 
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terband transition in the neighborhood of vaheys that 
exhibit circular dichroism. 

These criteria then lead us to manganese 
chalcogenophosphates, MnPXa, X = S, Se in 
monolayer form. Manganese chalcogenophosphates are 
layered crystalline materials, in which the interlayer cou- 
pling is the relatively weak van der Waals interactions. 
In principle, all van der Waals bonded layered com- 
pounds could be thinned to the monolayer limit by the 
micromechanical exfoliation technique. [15] |16] Hence, 
monolayer MnPXs in all likelihood can be produced. 
As shown in Fig. 4a, each unit cell in the monolayer 
MnPXa is composed of two Mn^+ and one [P2X6]^~ 
cluster, the latter isostructural and isoelectronic to 
molecular ethane (C2H6) in the staggered conformation. 
The most crucial feature of these compounds is that 
each Mn^+ ions, assuming an 5 = | high-spin state, is 
antiferromagnetically coupled to its three NN's, forming 
the afm honeycomb lattice. 
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FIG. 4: Structure and antiferromagnetism of monolayer 
MnPXs. (a) structure showing the unit cell. Purple spheres 
are Mn, yellow X, and gray P. Some of the computed bond 
lengths are: P-P = 2.22 A, P-S = 2.04 A for MnPSs, and P-P 
=2.24 A, P-Se = 2.22 A for MnPSes. (b) Spin densities in 
one unit cell, presenting the antiferromagnetic configuration 
(the isosurface of 0.4 e/A^). Left: up spin; right: down spin. 



Then we use density functional theory [17 calculations 
within the generalized gradient approximation (GGA) 
[m [19] to investigate the basic electronic structure of 
monolayer MnPXa. The projector- augmented wave po- 
tentials are used, as implemented in the Vienna Ab initio 
simulation package. [20l [21] A plane wave cutoff of 600 
eV and a Monkhorst-Pack k-point mesh of 30 x 30 x 1 per 
reciprocal unit cell are adopted. Vacuum slabs at least 15 
A thick are inserted between monolayer MnPXa to min- 
imized interaction between periodic images. Structure 
optimizations are performed with a convergence thresh- 
old of 0.01 eV/A on the interatomic forces. To account 
for the magnetic structure of the divalent transition metal 
Mn, we use a GGA+ U approach to describe the on-site 
electron-electron Coulomb repulsion. [22] The value of the 
isotropic U is set to 5 eV, as suggested by a previous as- 
sessment of this parameter for the divalent Mn^+. [23] 



As a calibration for the choice of the band gap of 2.4 
eV is obtained for bulk MnPSs, to be compared with the 
experimental band gap of bulk MnPSs of 2.7 eV. [24 For 
both bulk and monolayer MnPXs, the antiferromagnetic 
order of MnPXs is indeed the more stable in our calcula- 
tions, compared to the no n- magnetic and ferromagnetic 
states by at least 45 meV per a unit cell. The magnetic 
moment on each Mn is computed to be about 4.6 Bohr 
magnetons, in good agreement with the experiments. [14] 
Therefore, the choice of U is reasonable and used in all 
calculations. We compute the interband transition ma- 
trix elements using density functional theory at the linear 
response level. [3 ^ 

The spin density isosurfaces for MnPSs are displayed 
in Fig. 4b. The spin densities for the up and down spins 
show little difference on the [P2X6]^~ framework. But 
the densities for two spins on Mn are indeed well sepa- 
rated and localized on the two Mn ions, providing the ex- 
pected spin-contrasting asymmetric potential (Fig. lb). 
The band structures of monolayer MnPXs are shown in 
Fig. 5a. We observe that the band widths near the band 
gap are quite narrow for MnPXs, between 0.17 and 0.63 
eV. The ratio of U to the band widths (a measure of 
U/t) for MnPXs is greater than 7, indicating that the 
afm insulating state is expected to be stable in the mono- 
layer limit. [12 MnPSs and MnPSes both show direct 
band gaps at K± (2.53 and 1.87 eV, respectively), which 
fall well within the optical range. This is advantageous 
for the realization of optical polarization of charge car- 
riers. Within the collinear treatment of magnetism, the 
Bloch states of two spins are degenerate everywhere in 
the momentum space. When treating the spins in the 
non-collinear formalism to account for possible spin-orbit 
(or spin- valley) interactions, the spins remain degener- 
ate, but band gaps become renormalized as expected. 
For MnPSs, the gap difference between Kj- is negligibly 
small. But for MnPSes the gaps at K± differ by up to 
~ 43 meV, offering a window for realizing the SHE with 
equilibrium n- or p-doping (c./. Fig. 2b). 

To assess the optical selectivity of valleys by circularly 
polarized light, we compute the spin-dependent degree 
of circular polarization, ?7^*^(k),[6l |7] between the top 
valence bands and bottom of conduction bands. This 
quantity is computed using the density-functional linear 
response approach, defined as. 



p(.)(k)|2_|p^;(k)| 



is) I 



P|^^(k)P + |Pi^'^(k) 



(5) 



where P^*^(k) are the interband matrix elements of, re- 
spectively, left- and right-polarized radiation fields for 
spin s at k, defined for a vertical transition from band 
n to band n\ as P^^^(k;n,n') = (n'^k^slpx ±ipy|nk, 5), 
assuming spin flip is absent in the optical processes. The 
value of 77^^^ (k) quantifies the relative absorption rates 
of left- and right-handed photon. As shown in Fig. 5b, 
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monolayer MnPXa is computed to have perfect circular 
dichroism at K±. For one spin component they show the 
valley- selective circular dichroism as previously found for 
M0S2, [Z] and for opposite spins they show a selectiv- 
ity r]^^\r) = 25r, entirely concordant with our theoreti- 
cal model. The computed selectivity decays as we move 
away from the high-symmetry K±, but yet substantial 
polarization is achievable owing to the sizable regions of 
non- vanishing selectivity. 




FIG. 5: Electronic structure of manganese chalcogenophos- 
phates from DFT calculations, (a) The band structures of 
MnPXa (X = S, Se) near the band gaps, (b) The momentum- 
resolved degrees of circular polarization of MnPXs (X = S, 
Se), ?7*^^'*(k), between to the top of valence bands and the 
bottom of conduction bands. Only the values of one spin is 
presented, as in our calculations the other spin takes values 
equal in magnitude but with opposite signs over the Brillouin 
zone. At F, the computed optical selectivity is non-zero. This 
is a numerical artifact because of indeterminacy in 77^^-* (k) in 
the presence of degeneracy (apart from the spin degeneracy) 
in the valance bands (see (a)). 

The optical selectivity of valley interband transitions 
has been discussed in earlier papers (see, e.g., refs. [6]- 
[8]). Here, we phrase it in the language of group theory. 
Central to the selectivity is the rotational symmetry, Cs, 
and the the associated pure rotational group, C3. The 
point group of a staggered honeycomb lattice (A and B 
sites are inequivalent) C^y. The symmetry at K± is, 
however, C3; i.e., this is an Abelian group that does not 
allow degeneracy (apart from spin degeneracy). As band 
degeneracy is absent, then when we operate Cs on an 
eigenstate at a k that respects Cs rotational symmetry, 
we have Caink) = exp(i(/?^k)|^k), where (/?nk = ^nk27r/3. 



The azimuthal quantum number, /^k, is an integer deter- 
mined up to modulo 3. Consider the following transform, 

(6) 

Clearly, the optical selection rule is modulo(/nk — 
^n'k, 3) = ±1 for left- and right-polarized radiation fields, 
respectively. This is a mere re-statement of the conserva- 
tion of angular momentum in the absorption of a single 
photon absent spin fiip. The azimuthal quantum num- 
ber of the Bloch state, /nk = + ^m, which is the sum 
of two terms: (1) II arising from lattice phase winding 
of the Bloch phase factor, e*^ '^, [6 and (2) Im = 
where mi is the magnetic quantum number of the local 
atomic basis (or, Wannier-like functions) that contribute 
an additional phase under Cs. The effect of the local 
orbital symmetry and orbital ordering at the two valleys 
has an important role in determining the selection rules, 
as demonstrated in previous work. [7 

We would also like to remark on how spin-valley op- 
tical selection rule develops as a consequence of sym- 
metry breaking in a Fermionic honeycomb lattice. The 
three-fold rotational invariance of the Hamiltonian im- 
plies P^\'K) = for either left- or right-handed (±) po- 
larization field, when conduction and valence bands are 
separately non-degenerate (this is not the spin degener- 
acy, and is ensured when the point group is Abelian). [7] 
The antiferromagnetic Hamiltonian of Eq. (1) violates 
both time-reversal symmetry and parity, but accommo- 
dates the joint operation, O = TV] that is, [H''^^^ O] = 
0. When operating on a Bloch state, O preserves its 
wavevector but inverts spin. Consequently, the Bloch 
state of the same k and opposite spin are degenerate. 
Moreover, O changes the chirality of light, which implies 
P^\Kr) = P^~'\Kr), to within a phase factor. This 
is concordant with the notion that the valley-spin index 
s • r constitutes a degree of freedom, which defines chiro- 
optical selectivity (as well as orbital moments). 

Of note, the proposed optical selectivity may be a con- 
venient assay for antiferromagnetic order in an extremely 
thin sample, e.g., a monolayer or few-layer sample. This 
is especially relevant to the recent interest in the yet-to- 
be-uncovered spin liquid phase on a honeycomb lattice in 
the intermediate U regime. [12 The neutron scattering 
technique commonly used to determine the antiferromag- 
netic order will become impractical for thin samples. The 
circularly polarized photoluminescence, on the contrary, 
works well for these atomically thin samples with suffi- 
cient sensitivity. [71 m [To] In the spin liquid phase, ground- 
state ffuctuations are incessant at T = OK, whereby 
obliterating the afm order and, hence, the optical selec- 
tivity will be absent. Entrance into the symmetry-broken 
afm ordered regime, on the other hand, is accompanied 
by the proposed optical selectivity. 

A non-trivial variation of the proposed Hamiltonian of 
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Eq. (1) occurs when the spin-orbit interaction becomes 
overwhelming, compared to the afm order, m. Clearly, 
when \5\ = m (see Fig. Ic), the band gap closes in one of 
the valleys. When |^| > m, band inversion of one valley 
will drive the system into a Z2 topological quantum spin 
Hall state. [26l|27] In the limit m = 0, the model reduces 
exactly to that of Kane and Mele [26 , which affords time- 
reversal symmetry-protected edge states. However, when 
m is finite but smaller than time-reversal is lost. The 
edge states of a ribbon sample will be connected by the 
O-symmetry. This topic and the strategy for tuning the 
ratio \5\/m in real materials are issues worth further pur- 
suit. 

We also note that bilayer or few-layer MnPXa can 
also have similar optical selectivity, similar to sponta- 
neous symmetry-broken few-layer graphene system. [28 
Different kinds of stacking and magnetic order of bi- 
layer MnPXa (X = S, Se) with intra-layer antiferro- and 
ferro-magnetism are evaluated in our calculations, where 
van der waals corrections within GGA are included by 
optB86b-vdwDF method. [H [30] Two AA-stacked bilay- 
ers with intra-layer antiferromagnetism are nearly de- 
generate in energy, and more stable than other stacking 
and magnetic orders. They differ by the inter-layer mag- 
netic configuration; that is, the interlayer couplings are, 
respectively, ferromagnetic and antiferromagnetic. By 
computing the momentum-resolved spin dependent de- 
gree of circular polarization, we find that the s • r-CD is 
present in ferromagnetically coupled bilayer, which has 
no inversion center in the magnetic space group. The 
chiral optical selection rule is absent in the case of antifer- 
romagnetic interlayer coupling, as expected. As the per- 
fect spin-dependent optical selectivity of both monolayer 
and ferromagnetically coupled bilayer MnPXa, which at- 
tests to the emergent electronic degree of freedom pre- 
sented in the foregoing model, it is therefore of inter- 
est to experimentally study these materials in monolayer 
and few-layer forms, to characterize their optoelectronic 
and transport behaviors and to explore its potential for 
application in novel operating paradigms for advanced 
electronics. 
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